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Abstract. A new method of computing the homomorphism induced by a continuous map 
in homology presented in [1] and [2] relies on computing coboundaries of cycles. This 
paper is devoted to a precise geometric construction of a coboundary of a given cycle in a 
prescribed rectangle and a description of the associated algorithm. 


1. Introduction 


In a recent paper [2], anew method of computing the homomorphism induced by a con- 
tinuous map f in homology is presented. The method is based on finding a rectangular- 
valued point-to-set map F whose graph contains the graph of f and then finding a chain 
selector of F. In this approach, the computability of the homomorphism in homology 
depends on the ability of solving the following problem: 

Given a q-dimensional cycle z supported in a rectangular set A, construct a (q + 1)- 
dimensional chain c, also supported in A, such that 


dc =z. (1.1) 


Since any rectangular set A is acyclic, a solution c always exists. We call it a cobound- 
ary of z and denote it by COB(z). Note that c is not unique in general. Finding c by direct 
matrix algebra methods results in a large number of algebraic operations, especially due 
to the fact that the involved matrices are not invertible and that computation must be done 
for integer coefficients, see [8]. Consequently, linear algebra leads to a high complexity 
of the final algorithm computing the homology of a map. 


* The second author was supported by NSERC of Canada. 
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We present, in this paper, a geometric construction which takes advantage of a rectan- 
gular form of the acyclic set A. The algorithm based on this method computes COB(z) in 
linear time with respect to the length of z, that is, with respect to the number of generators 
in the carrier of z. We briefly mentioned in [2] how this construction can be done but, 
considering its importance, it seemed to us relevant to give a precise description of the 
algorithm and a proof of the correctness of the formula for COB(z). 

Our main motivation for studying (1.1) comes from dynamical systems where the 
homomorphism induced in homology by a continuous map carries important information 
about the dynamics of a map. We refer the reader to [2], [9], and references therein for 
further discussion on that topic. 

Another potential motivation for studying this problem is in the fact that it is related 
to the problem of finding minimal surfaces of closed curves (see [3]) and, in particular, 
Seifert surfaces (see Chapters 4 and 5 of [13]). This is a hard and beautiful problem 
intensely studied by methods of differential geometry, calculus of variations, and com- 
binatorial geometry, and we do not claim that our construction will make a miracle here. 
In order to adapt our construction for the study of minimal surfaces, one must generalize 
it for simplicial triangulations since cubical surfaces will always be far from minimal. 
For example, the shortest path joining two opposite vertices of a unitary square [0, 1]? 
is the diagonal of length /2 whereas our algorithm in the unitary cubical grid would 
give a sum of two edges of total length 2. Generalizing our construction to simplicial 
triangulations is a work in progress. If this is done, the two problems can be related as 
follows: 

Let C be a polygonal closed curve with consecutive vertices 


vo, U1, v2, ae Un-1s Un = vo 


which approximates a given smooth curve K. The first rough attempt would be to take a 
convex hull A (it is an acyclic set) of all vertices of C and triangulate it. Then the support 
of the cycle 


z := [vo, vi] + [v1, v2] +--+ + [Up-1, vo] 


of the simplicial complex of A is C. If we solve (1.1), the support S of c will be a two- 
dimensional surface whose one-dimensional skeleton contains C. Obviously we have no 
guarantee that our surface is close to a minimal surface of K but maybe the search for 
an optimal triangulation and the best solution of (1.1) can be somewhat refined. 

To the authors’ best knowledge, the problem (1.1) in arbitrary dimension has not been 
studied from the computational point of view prior to [1] and [2]. For recent results on 
computing homology of spaces, we refer the reader to [4], [5], [10], [12] and [18]. For 
an overview of problems of computational topology, refer to [6]. 

Since the first submission of this paper, important progress has been made in the 
direction of designing programs computing the homology of spaces and maps. A software 
based on the algorithms presented in [10] has been developed by Pilarczyk [19]. A parallel 
work based on [12] has been done by Kalies and Watson. An algorithm for solving (1.1) 
by the above discussed linear algebra methods has been found by Mazur and Szybowski 
[15]. The authors of that paper also wrote a computer program for the homology of 
a multivalued map F with two options to choose: our geometric algorithm or their 
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algebraic algorithm. Their experiments seem to confirm the efficiency of the geometric 
construction. That program is included in the Pilarczyk folder [19]. 

The paper is organized as follows. In Section 2 we briefly summarize preliminary 
definitions related to cubical grids and cubical chain complexes. A more systematic 
exposition can be found in [11]. In Section 3 we start from a simple example aimed at 
providing the reader with geometric intuitions behind the recursive formula for COB(z). 
We next state the formula and prove Theorem 3.2 which states that OCOB(z) = z and 
thus it affirms the correctness of the construction. Although the construction has a clear 
geometric flavor, the proof of the Theorem 3.2 is not trivial and it constitutes the main 
achievement of this paper. In Section 4 we give some more examples illustrating the 
construction. In Section 5 we describe the algorithm and discuss its complexity. 


2. Basic Concepts 


The notions of grids, representable sets, and representable maps are presented in their 
general formulation in [16]. Representable sets and maps for convex polyhedral grids are 
discussed in [1] and [2]. We recall in this paper the most significant case of cubic grids in 
R”. The cubic grid Ex of mesh 1/k, k = 0, 1,2,..., is the collection of all cubes e € R” 
of the form e = (1/k)(x + I, x h x ... x In) where x € Z” and J; is either the open 
interval ]0, 1[ or the singleton {0}. We call the elements of E elementary cubes. Note that 
R” is the disjoint union of e € Ex. There are many advantages of such decompositions 
over simplicial triangulations, for example, the product of two elementary cubes is an 
elementary cube which is not true about simplexes. A subset X of R” is said to be 
representable over E; or, shortly, k-representable, if there exists a finite €’ C Ep such 
that X = |) €’ and X = X. In particular, X is a finite polyhedron. The set £’ is uniquely 
defined and denoted by &(X). The family of all k-representable polyhedra in R” is 
denoted by Rg. 

Any cubic grid E% can be brought to the unitary cubic grid € = € by simple rescaling 
x œ> kx, x € R”, which is an isomorphism of R”. This isomorphism maps the set Rg 
to R = Rı. Ina similar way, multiplying by different factors on each coordinate axis 
permits rescaling general rectangular grids to get €. Grid refinements are only needed 
for purposes of numerical approximation of sets and maps. In this paper we always work 
in one fixed grid, so it is enough to present the definitions and results for the unitary 
cubic grid E. 

Given X € R, the cubical chain complex of X is the pair of sequences 


C(X) = (Cg (X) qez, (8g: Ca (X) > Cq-1(X) qez) 


defined as follows. Each C; (X) is a free abelian group generated by the set of elementary 
cubes €4(X). We put C,(X) = 0 if g < 0 org > n. Thus the elements of this group, 
called chains, are formal linear combinations of elementary cubes. The carrier of a chain 
c = } ce; is the cubical set defined by 


lel = (J i: ci #0}, 


where g; is the closure of e;. 
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The map 0,: Cg(X) —> C4-1(X), called the boundary map, is a group homomorphism 
defined below on elementary cubes and extended by linearity to all chains. First, for the 
cell 74 = ]0, 1[% C R4, we put 


Ail? = (1, ..., Xi—1, 0, Xi+1, - -< Xq): x E€ 11} 


and 
BI? = {(x1,.--,Xj-1, 1, X41, ven Xg)! XE I}, 


The faces A;/% and B;I‘ are called, respectively, the front i-face and the back i-face 
of 11. We let 


q 
əl = X (Dİ (Ail? — Bi’). 
i=1 


Now let q-cell e € €7(X), X C R” be of the forme = x +1], x h x --- x I, where 
x € Z” and J; is either the open interval ]0, 1[ or the singleton {0}. The formula for d,e 
is more complicated than that for 77 since some of the intervals Z; may be singletons 
and, in this case, the jth coordinate does not contribute to the sign alternation. It is easily 
seen that e is the image of J? under the affine-linear map T given by the formula 


Y1, Y2,+++5 Yn ee . 
7 xj + Yj if I; =]0, 1[. 

The map T is composed of an insertion of zeros on chosen (n — q) coordinates 
corresponding to the indices j such that J; = {0} with the translation of the origin to the 
vertex x € Z”. We put Aje = T(A;I‘%), Bie = T(B;I‘), and define 


q 
age = X- (—1)' (Aje — Bie). (2.1) 


j=l 


Note that ð = 0 since C-ı(X) = 0. For simplicity of notation we avoid indices and 
write ð for 0, whenever the index is clear from the context or if the statement is true for 
all g. The most important property of the boundary map is that 


a03=0, (2.2) 


i.e. 0g—1 © ðq = 0 for all q. For the proof that our cubical boundary map has this property 
we refer to [11] and [14]. A somewhat similar construction is also given in [7]. 

A chain z € C,(X) is called a cycle or, more precisely, a q-cycle if dz = 0. A chain 
z € C,(X) is called the boundary if there exists c € Cq+1(X) such that dc = z. The set 
of all g-cycles is the subgroup Z, := ker 0, of C,(X), while the set of all boundaries is 
the subgroup B, := im 0,,;. Equation (2.2) implies that im 0,4; C ker 0, and hence the 
quotient group H,(X) := ker d,/im 0441, called the qth homology group of X, is well 
defined. The homology of X is the sequence H,,(X) := {H,(X)}. 

For a reader who is not familiar with homology, we should mention much information 
can be extracted just from the dimension f, of the free component of the group H,(X). 
This number is called the qth Betti number of X. Thus fọ is the number of connected 
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components of X, 6; is the number of nontrivial fundamental (i.e. generic) loops in 
X, and, if X C R?, B> is the number of “holes” in X (think about holes in a brick of 
Ementaler cheese). 

A representable set X is called acyclic if it has the same homology as a point. The 
homology of a point p, however, is not all trivial: 


if 
minafe fate 


In order to associate trivial homologies with acyclic spaces one introduces reduced 
homology H,(X) as follows: H,(X) = H,(X) if q # 0 and H,(X) = kere/im ð}, 
where the map e: Co(X) — Z, called the augmentation map, is given on vertices 
(zero-dimensional elementary cubes) by ev := 1 for all v € Eo(X). 

The reduced q-cycles in Zs (X) are the same as cycles if q # 0 and they are elements 
of Zo(X) = kere if q = 0. It is easily seen that the reduced 0-cycles are generated 
by chains v — u where u, v are vertices. If X is acyclic, then all reduced cycles are 
boundaries. That implies, in particular, that any two vertices of X can be connected by 
a path of edges. For more details see, e.g. [17]. 

We finish this section with the remark that any representable cubical set X can be 
triangulated and its cubical homology defined above is isomorphic to the classical sim- 
plicial homology. Conversely, a polyhedron P is not necessarily a cubical set but it is 
homeomorphic to one, thus the two homology theories are equivalent. For details, refer 
to [11]. 


3. Coboundary Construction 


Let then A C R” be a representable rectangle, i.e. a Cartesian product of intervals with 
integer coordinates. Given a reduced cycle z € Zs (A), we construct a chain c € Cy4)(A) 
such that 0,41¢ = z. We denote this procedure by c := COB(z) where COB stands for 
the word coboundary. We construct in a general manner an operator 


COB: Cy (A) =— Capir); 


by defining COB (o) for allo € C4 (A). 
We start from a simple example. 


Example 3.1. Consider the clockwise oriented cycle z = eg + €7 + e6 + e5 — e4 — e3 — 
ez — e1, in R?, where each e; is a positively oriented interval as indicated in Fig. 3.1. 
Then 


z€C,(R(z)) where R(z) = [0, 2}. 


Let xı be the canonical projection onto the x2-axis and let R,(z) = {0} x [0,2] be 
the image of R(z) under this projection. For each interval e; that is not projected to 
itself or a point, define [77 (e;), e; ] to be the formal sum of the unit squares in which e; is 
projected down to 77 (e;), otherwise define [77 (e;), ei] = 0. It follows from the definition 
that [771 (e1), e1] = [m (e2), e2] = [m1 (7), e7] = [71(@s),e5] = [71 (es), eg] = 0, 
[7 (e6), e6] = E4, [711 (€4), e4] = E3 + E4, and [7]; (e3), e3] = Ey + Ep. 
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Fig. 3.1. Coboundary construction. 


Now define a 2-chain COB(z) by replacing each e; in the formula of z by [7 (e;), ei] 
and keeping the same coefficients, i.e. 


COB(z) = 0 + 0 + E4 + 0 — E4 — E; — E; E = — (E; + E2 + E3). 


It is clearly seen that z is the boundary of COB(z). Here, if one takes into account the 
orientation, then the projection of z onto the x2-axis is 0. One asks the question whether 
the definition above is still correct when the projection is not trivial. Although the answer 
is negative, the formula above is still the main ingredient for the correct construction we 
provide in the rest of this section. 


Let o € C,(A). One can write o = eee ciei where e; € E1(A) and c; Æ 0, Wi. 
Each e; can be expressed in the form e; = ((ai)1, (bi)1) X ++- X ((ai)n, (bi)n), where 


_ fa, of if aÆb, 
ab = fia if a=b. 
The endpoints of intervals (a;);, (b;), are integers, either (b;), = (aj), + 1 or (bj), = 
(a;),, and the first equality holds exactly for q values of k. Let then R (o) be the smallest 
representable rectangle in A containing the carrier |o| of ø. It is easily seen that 


R(o) = [m1, Mı] x [m2, M2] x «++ x [Mmn, Ma), (3.1) 


where mg = min; {(a;),} and Mg = max; {(b;),}. The dimension of R (o) is the number 
of nontrivial intervals in formula (3.1). 

Let kı be the first nontrivial coordinate of R(o), i.e. the smallest integer k with the 
property mg + M, (such an integer exists unless the chain ø is trivial). Consider now 
the rectangle 


Ry, (0) = {x € RG): Xk = me} 
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contained in the hyperplane orthogonal to the first nontrivial direction kı of R(o). Let 
Tk: R(o) —> Rx, (0) be the canonical projection defined by 


Xj if ixk, 


Mk, otherwise, (3.2) 


{1,0}; = | 


and let z4: C,(R(o)) —> C4(Rp (0 )) be the chain homomorphism induced by 7g, . The 
chain homomorphism sry is defined at level q on any g-elementary cube e = (a1, b1) x 
+++ X (an, bn) € Ca (R (0 )) by the formula 


0 if by, = ax, +1, 
Tyle) = (ay, bı) Kr X (ax,-1, Dk, -1) x {mk} X X (an, bn) 
if De, = k. (3.3) 


We construct, by induction with respect to the dimension d of R (o), a (q + 1)-chain 
COB(o). Obviously the dimension of R(øo) is at least q so the induction starts from 
d=q. 

For d = q, set COB(c) = 0. 

Let d > q. Suppose the construction is done for dimensions up to a certain d > q 
and now let dimR(o) = d + 1. Fore = (a), b1) X +--+ X (an, bn) in EY(R(G)), we define 


(0) if male) =0 or azg(e) =e, 
by, —1 Aeh 
[74(€),€] = } Doing, @1 b1) X0 X (ani bam) X (ED) X++ X (an, bn) 
otherwise. (3.4) 


Note that, by definition, [7(e), e] is a (q + 1)-chain in A. If r(e) Æ 0 and e is not 
contained in Rx, (0), then 


I[7r#(e), ell = (ai, bı) Kr XK (ax,—-1, bx,-1) x (mx,, bki) Kr X (an, bn) (3.5) 


is the (q + 1)-dimensional rectangle through which e is projected down to zrg(e). 

The (q + 1)-chain COB(z4(o)) is well defined by induction hypotheses since z(o) 
is a g-chain contained in Rg (o), which is a representable rectangle of dimension d. We 
then define 


N 
COB(a) := COB(r4(o)) + J cilza(e:), ei]. (3.6) 
i=l 
We may now state the main result of this paper. 


Theorem 3.2. Letz = Sr cie; be a reduced q-cycle whose carrier is contained in 
A. Then 


0g+1COB(z) = z. (3.7) 


Proof. We proceed by induction on dimension d of R(z). As in the construction of 
COB(z), we may start the induction from d = q. 
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Ifd =q. By definition, COB(z) = 0. On the other hand, we have Cori (RG) = = 0, 
then B,(R(z)) = Oand Z, (R(z)) = 0 since A; (R(z)) = 0. Knowing that z € Z g(R(Z)), 
it follows that z = 0, and ‘anally 0g41COB(z) = z. 


Let d > q. Suppose that for dimR(z) < d, (3.7) is satisfied. Now let dimR(z) = d+1. 
We begin by establishing the following result. 


Lemma 3.3. We have 


N 
Oq+1 (dre cjl#(e;), e:i] i) = z — mahz). 


i=l 


Proof of the Lemma. We proceed by induction with respect to the length M = Mp, — 

x, Of the first nontrivial side of the rectangle R (z) (which is the projection of R (z) onto 
the xg -axis). Note that the transition from M to M + 1 contains the arguments that prove 
the case M = 1. When the length is M + 1 = Mg, — m,,, the rectangle R (z) is made up 
of M + 1 slices, each of length 1, along the x,;,-axis. Consider 


R' = [m;, Mı] x [m2, M2] x -+ x [mp , Mk, — 1] x +++ x [mn, Mn], 


the rectangle made up of the first M slices of R (z). Define the retraction r: R(z) —> R’ 
by setting, for x € R(z), {r(x)}; = x; if i Ak, and 


XK, if xn SMa- 1, 
ræk = a M, -1 if xp, > My, -— 1. 


The chain homomorphism ry: C,(R(z)) —> C,(R’) induced byr at level q is defined 
on a q-elementary cube e = (a1, b1) X ++- X (an, bn) E Cq(R(z)) by the formula 


0 if akı = My, — 1, be, = Mk, 
ry(€)= į (a1, By) X ++ X (ak1, Bk 1) X (Mg, —1} X- --X (an, bn) if ag, = Mp = dy, 
e otherwise. (3.8) 


Set z’ = rx(z) = oG cirą(ei), the image of z by ry. By definition, z’ is a reduced 
cycle in R’. It is clear that for any x € R(z), mx, (r(x)) = mx, (x) (see formula (3.2)) and 
therefore Rg, (z) = Rx, (z’). One can also invoke formulas (3.3) and (3.8) to check that 
Ty(ra(e)) = ma(e) for any e € E4(R(z)). It follows that w4(z) = m#(z’). 

The carrier of the reduced cycle z — z’ in C,(A) is contained in the last slice R” of 
R(z) given by the formula 


R” = [m,, Mı] x [m2, M2] x «++ x [Mp — 1, Mi] x +++ x [mn, Mn]. 
At this point, we decompose each (q + 1)-chain [7tx(e;), e;] as 
[#(e;), ei] = [r#(e;), re(e:)] + [re(ei), ei], (3.9) 


where [zrx(e;), r#(e;)] is defined as in formula (3.4) by noticing that zrx(e;) = m(r#(e;)) 
and 
0 if rg(e)=O0O or rg(e) =e, 
[r#(ei), ei] = § (a1, b1) X +++ X (Gy-1, Dg,-1) X (Mg, — 1, Mk) X +++ X (an, bn) 
otherwise. (3.10) 
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The carrier of the (q + 1)-chain T = pe ci[rx(e;), ei] is contained in R” and one 
can write 


Me 
oe 


cilmg(e;), ei] = cilg(ei), raele] +T 


i=l 


il 
z 


= $ cilrglrglei)), ree] +T. (3.11) 


i=l 
We first prove that 
Ogi =z- 7. (3.12) 


This also proves the case M = 1. Note that if M = 1, we have zy = rg, R(z) = Rp (Z), 
z’ = T(z), andl = as cj[7#(e;), ei]. Therefore, proving formula (3.12) will prove 
in case M = 1 that 


N 
Og+1 (Ea cilg(ei), e:i] i) = Z — mahz). 


i=1 
To prove formula (3.12), it is sufficient to establish that for any cell e € E4(R”), the 
coefficient of e in z — z’ is equal to its coefficient in 0,411. With this in mind, we define 
on C,(A) x C,(A) a bilinear form by setting, for any two cells e, o € E1 (A), 


ees 1 if o=e, 
met 10 otherwise. 


It is then enough to verify that 
(e, Ig4iT) = (e,z — 2’), Ve € E1(R"). (3.13) 


Let Px, be the projection of R (z) onto the x;,-axis. We can then consider three types 
of q-cubes in R”: 


1. The set A of q-cubes A such that P, (A) = {Mp — 1}. 
2. The set 6 of q-cubes B such that P, (B) = {M;,}. 
3. The set C of q-cubes C such that Py, (C) = (Mk, — 1, Mp). 


The (q + 1)-cubes in € 4+1 (R") are in one-to-one correspondence with the elements 
of A and with the elements of B. Indeed, for each 


Q = (a, b1) x ++- x (My, — 1, Mk) x +++ xX (an, bn), 
we can find a corresponding 
B = (aj, b1) X +++ X {Mpk} X +++ X (an, bn) E B 
and 
A = (aj, b1) X +++ X {Mk, — 1} X -+ X (an, bn) = re(B) € A. 


One can see that, for any B € B, the q-cube A = ry(B) € A. More precisely, the 
restriction of rz to B is a bijection from $ onto A. For a given B € B, we denote by 
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[A, B] the (q + 1)-cube obtained by projecting B down to A. It is easily noted that if 
Q € €4*'(R”) and B € B such that Q Æ [A, B], then 


(B, dQ) = (A, dQ) =0. 


Note that, since B is a front face of [A, B] and A a back face of it in the sense given in 
the definition of 0, it follows from the choice of kı that 


(B, o[A, B]) = 1, (A, o[A, B]) = —1 (3.14) 


for any any B € B such that (B, z) 4 0. 
We deduce from what precedes that for each B € B, one has 


(B, ei) = ([A, B], [rg(e;), ei]), vi. 


Thus, 
(B, z) = ([A, B], T) = (A, z^ — (A, z), YB €B. (3.15) 


The latest equality is due to the fact that A could be present in the formula of z. Thus its 
coefficient in z’ is the sum of its coefficient in z and the coefficient of B in z (see Fig.3.1 
where the cell eg = r4(e6) = r (e4) and therefore its coefficient in z’ is 0). 

Then we show the identity (3.13) for elements of A, 6, and C: 


1. Let B € Band b = (B,z — z’). Since B does not intersect R’, it follows that 
(B, z') = 0. Thus b = (B, z}. One can then deduce that 


(B, OV) = b(B, o[ryx(B), B]) = b. 


2. Nowlet A € Aand B € B be such that ry(B) = A. Leta = (A, z) andb = (B, z). 
By formula (3.15), we have 


b= ([A, B], r) = (A, z) e (A, z). 
It follows that (A, z — z’) = —b. On the other hand, 
(A, ðI) = b(A, Ogi lA, B1) = —b. 


3. Finally, letC € Canda = (C, z—z'). By definition (see formula (3.8)), r#(C) = 0. 
Therefore, (C, z’) = 0, hence œ = (C, z}. Since l = yS ci[lrg(ei), ei], it follows 
that 


N 


(C, dg4iT) = (c.¥ Dek Ogi (H(i), ei] = Yo cilC, Ig41 (Irale), e;)). 


i=l 


Obviously, if the index i is such that [rg(e;), ei] = Oor C is nota face of [rx(e;), ei], 
then (C, 0,+1[r#(e), eil) = 0. Then let J be the set consisting of the indices 
i € {0,1,..., N} for which [rs(e;), ei] # 0 and C is a face of [rg(e;), e;]. For 
i € J,denote bya; = (C, 0,41 [r#(e;), e;]) the coefficient of C in 0,41 ([r#(e;), e]). 
Then 


(C, 3T) = YC ciai. 


ieJ 
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Let now F be the (q — 1)-face of C located at the level x,, = M,,. It is easily seen 
that F is a common face of C and all e;,i € J. It follows that 


0 = (F, 3gðq+ı(lrz(ei), ei))) = (F, 3g (aC + e;)) 


= (F, ðei) +a; (F, jC), Wied. 


On the other hand, since z is a cycle and using the precedent equality, we have 


ieJ ieJ 


0 = (F, 394z) = J ci (F, dge;) +a (F, 34C) = (- = Yon] (F, 3C). 


Hence a = pe c;a;, which proves formula (3.13) for elements of C. 


Up to now we proved that 0,4;1F = z — z’, and 


N 
Og41 (> cilra(ei), 7) = Z — T(z), 


i=l 
for M = 1. Let M > 1 and suppose that the length of the first nontrivial side of R (z) (in 
the xg, direction) is M + 1. We recall from formula (3.11) that 


N N 
S cilrglei), ei] = J cilrglrgle)), ree +r. 
i=l i=l 
The smallest rectangle R(z’) containing the carrier of z’ and the carrier of the (q + 1)- 
chain Sh cilt#(re(e;)), r#(e;)] is contained in R’. The first nontrivial side of R(z’) is 
also in the xg, direction and is at most of length M. By the induction hypotheses, formula 
(3.12), and the fact that z4(z) = s4(z’), we can finally write 


N 
Og+1 (de Ci [7x (e;), ei] 1) 


1 


N 
ðq+1 (£ cilt (rg(e;)), reni) + q(T) 


i=1 


= z — mz) +z- 
= Z — T(z). 


This completes the proof of the lemma. 


The smallest rectangle R (7r4(z)) containing the carrier of 74(z) is contained in Rx, (z), 
the rectangle obtained by projecting R(z) onto the hyperplane orthogonal to the direction 
Xk. It follows that the dimension of R(sr(z)) is at most d. Then, from the induction 
hypotheses and Lemma 3.3, we deduce that 


N 
0g-+1(COB(z) ðg+1ı (COB (t4(z))) + Og-41 (> ci[tg(e;), ei i) 


i=1 


Ty(Z) +z — T(z) 


= Z, 


and the theorem is now proved. 
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} 


X 


Fig. 4.1. Coboundary of a 0-cycle in R?. 


4. Examples in R? 


The following examples provide more insight into the construction given in Sec- 
tion 3. 


Example 4.1. Consider in R? the reduced 0-cycle z = e2 — e1, where e; = {(0, 0, 0)} 
and e2 = {(1,1, 1)} (Fig. 4.1). Here R(z) = [0, 1]? is of dimension 3 and its first 
nontrivial direction is 1. 

Call xı the projection onto the x2x3-plane. The rectangle R(z) is projected onto 
Rı(z) = {0} x [0,1]? and z is taken to the cycle zı = (71)y(z) = e’2 — e1, where 
e’ = {(0, 1, 1)}. The dimension of R; (z) is 2 and its first nontrivial direction is 2. If 72 
is the projection onto the x3-axis, then R; (z) is projected onto the rectangle R2(z1) = 
{0} x {0} x [0, 1] of dimension 1 and with x3 as the first nontrivial direction. The cycle 
zı is projected onto the cycle z2 = (72)4(z1) = e”2 — e1, where e” = {(0, 0, 1)}. Finally 
the projection onto the origin yields the rectangle R3(z2) = {0} x {0} x {0} of dimension 
0 and z? is taken to the cycle z3 = e; — e; = 0. It follows from the construction that 
COB(z3) = 0 and therefore 


COB(z) = COB(z1) + [e’2, e2] 
COB(z2) + [e"2, e'2] + [e’2, e2] 
= COB(z3) + [e1, e”2] + [e"2, e'2] + [e"2, e2] 


[e1, e”2] + [e"2, e'2] + [e’2, e2]. 
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x 


4 


Fig. 4.2. Coboundary of a 1-cycle in R?. 
Example 4.2. Consider in R? the 1-cycle 
6 
z= X ciei =e; — 2 + e3 + €4 + €5 — e6 
i=1 


where each e; is a positively oriented interval as indicated in Fig. 4.2. 

The smallest rectangle that contains the carrier of z is R(z) = [0, 1]°. Its dimension 
is 3 and its first nontrivial direction is 1. The projections zr, and zz are defined as in the 
preceding example. Consider the open squares E;, i = 1, ..., 6, where E is the base of 
the cube R(z), E2 is the face having e; and e6 as sides, E} is the face having e2 and e3 as 
sides, and E6, E4, Es are the respective opposite faces to E1, E2, and E3. It follows that 
Ri (z) = {0} x [0, 1]?. Moreover, m1 4(e1) = e1, Mig(€2) = e2, Tig(e3) = T14(e€6) = 0, 
i4(€4) = eg, and m1 4(e5) = e7. Finally, m14(z) = —e; — e2 + eg + e7. It is easily seen 
that the projection of ;4(z) onto the x3-axis is W24(714(z)) = 0. On the other hand, 
since 


[7r24(€1), €1] = [T24(e8), eg] = [724(e2), e2] = 0 
and [7r24(e7), e7] = E; we have COB(14(z)) = E1. Finally 


6 
COB(z) = E, + X cibrig(ei), ei] 
i=l 


= E£,-0-0+0+ £4+ Es; —0 
= Ei + E, + Es. 


Remark 4.3. Obviously, there are many chains c such that dc is a given cycle. We may 
obtain alternative constructions by projecting z onto different faces of the cube R(z), 
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but then the signs may change in (3.14), thus care must be taken about the alternation of 
sign in the formula for a boundary of a cube [2], [11]. 


5. The Algorithm 


We now discuss the formal algorithm permitting the construction of a coboundary 
COB(z) for a reduced cycle z, namely, a chain COB(z) such that ODCOB(z) = z. 

Let z = ue ciei be a qg-reduced cycle in R” where 0 < q < n and each particular 
ei is a (generator) cube of dimension q in R” that is written in the form 


ei = (4), (bi)1) X +++ X Cas Bd a); 


where only q intervals among {((a;);, (bi) ;) sai are nondegenerate. 


The cubes will be represented by the type CUBE. The dimension of an element of type 
CUBE is internal to its representation. The type CHAIN will be used to represent a list of 
elements of type CUBE of the same dimension together with their coefficients (integers). 
The dimension DIM (integer) of an element of type CHAIN is the common dimension of 
all its elements of type CUBE. A chain is also characterized by its length LENGTH which 
is an integer giving the number of cubes in the chain. We do not go into the details of the 
definition of each of these representations here. The smallest rectangle containing |z| is 
given by 


R(z) = [m;, Mı] x [m2, Mp] Xr X [mn, Mal, 


where mg = min; {(a;);} and M, = max;{(b;);,}. The dimension of R(z) is defined to 
be the number of indices i € {1,..., n} for which m; ~ M; and is less than or equal 
to the dimension of the ambient space R”. The first nontrivial direction of R(z) is the 
first index i for which m;  M;. We therefore use the type RECTANGLE with associated 
dimension DIM and first nontrivial direction FIRST (necessary to define the projection) 
in its representation. 

We can create a rectangle from a given chain z, 


RECTANGLE(z: CHAIN); 


in at most nN steps where n is the dimension of the ambient space R” and N is the 
number of cubes in the chain z, i.e. the length of z. 
We make use of two functions, namely, 


PROJECTION(R: RECTANGLE; e;: CUBE): CUBE; 


and 


LIFT(R: RECTANGLE; e;: CUBE): CHAIN; 


that are simply defined according to the definitions in formulas (3.3) and (3.5) respec- 
tively. These two functions when called run in a constant time provided the output of the 
function RECTANGLE is presented in an appropriate way. For more details on algorithms 
see [20]. 
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Algorithm 5.1 


function COB(z: CHAIN): CHAIN; 


begin 
R := RECTANGLE (z); 
d:= DIM(R); 


q := DIM(z); 
N := LENGTH (z); 


if d=q then 
return 0; 
else 


for i:=1 to N do 
Oi := PROJECTION (R, ei); 
Ei = LIFT (R, 0;, €i) ; 
end for; 
o := D CiOi} 
E= SG ci Ei; 
return ŒE + COB(0); 
end if; 
end; 


The recursive algorithm runs d — q + | times, and a single call costs O(N) where 
N is the length of the cycle z, q is the dimension of z, and d is the dimension of the 
rectangle R(z). Since we work in a fixed dimension n, then d — q + 1 < n + 1 and the 
algorithm runs in linear time in the number of cubes in the initial chain. 


Remark 5.2. We mentioned in the Introduction, that it is possible to solve the equation 
dc = z by row reduction of noninvertible matrices. When this approach is chosen, the 
complexity of the algorithm is about O (N 3) in case of field coefficients and about O (N5) 
when integer coefficients are assumed. 


Remark 5.3. Our construction of c = COB(z) is not minimal in the sense of the 
number of generators in the support of c. That could be improved by a judicious choice 
of projection faces in each step. This however will substantially increase the complexity 
of the algorithm. 
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